Molecular dynamics simulation study of the high frequency sound waves 
in the fragile glass former ortho-terphenyl. 
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Using a realistic flexible molecule model of the fragile glass former orthoterphenyl, we calculate via 
molecular dynamics simulation the collective dynamic structure factor S(Q,u), recently measured in 
this system by Inelastic X-ray Scattering. The comparison of the simulated and measured dynamic 
structure factor, and the study of the S(Q,u) in an extended momentum (Q), frequency (w) and 
temperature range allows: i) to conclude that the utilized molecular model gives rise to S(Q,u>) 
in agreement with the experimental data, for those thermodynamic states and Q values where the 
latter are available; ii) to confirm the existence of a slope discontinuity on the T dependence of 
the sound velocity that, at finite Q's, takes place at a temperature T x higher than the calorimetric 
glass transition temperature T g ; ni) to find that the values of T x is Q-dependent and that its Q — »0 
limit is consistent with T g . The latter finding is interpreted within the framework of the current 
description of the dynamics of supercooled liquids in terms of exploration of the potential energy 
landscape. 

PACS numbers: 64.70.Pf, 71.15.Pd, 61.25.Em, 61.20.-p 



I. INTRODUCTION 

The understanding of the dynamics of the supercooled 
liquids flW and of its relation with the glass transition 
has received great attention in the last few years. Dif- 
ferent details of the dynamics of supercooled liquids and 
glasses seem to be on the way to be settled, while other 
are still obscure. As an example, considering the dynam- 
ics of the density fluctuations which are the main issue 
for this paper, on one side the connection between the 
glass transition phenomenology and the long time dy- 
namics (structural rearrangement) has been almost clari- 
fied, on the other side the effect of the structural arrest on 
the high frequency collective vibrational motion (sound 
waves) is much less clear. 

From the experimental point of view the study of the 
time behavior of the density fluctuations is made via the 
determination of the dynamic structure factor S(Q,uj), 
i. e. the power spectrum of the Q component of the 
number density PQ(t). In the hydrodynamic regime this 
quantity is measurable by laser light scattering, a well 
established technique (Q «0.03 nm _1 corresponding to 
frequency in the GHz range). Thanks to this technique 
much is known about the collective dynamics in the 
nanosecond time scale. The limit of high frequencies (pi- 
cosecond time scales) is traditionally the realm of the 
inelastic neutron scattering technique. This technique, 
however, suffers of strong kinematics limitations, that 
prevent the possibility to investigate the mesoscopic Q 
region (Q « 1 v 10 inn -1 ) that is particularly interest- 
ing as it is the region where the transition from a gen- 



uine collective behaviour of the density fluctuation fades 
in the single particle dynamics. Such kinematics limita- 
tions have been recently overcome by the Inelastic X-ray 
Scattering (IXS) technique. Using this techniques it has 
been possible to firmly establish few common features 
of the collective dynamics in glasses in the mesoscopic 
region In particular, beside from specific quantita- 
tive difference among different systems, all the investi- 
gated glasses show some qualitative common features that 
can be summarized as follows: i) there exist propagating 
acoustic-like excitations up to a max Q value (Q m ) given 
by Qm a ~ 1 -j- 3 showing up as more or less defined Bril- 
louin peaks at fl(Q) in the S(Q,u>) (here a is the average 
inter-particle distance). The specific value of Q m a re- 
sult to be correlated with the fragility of the glass ; ii) the 
slope of the (almost) linear £l(Q) vs Q dispersion relation 
in the Q — > limit extrapolates to the macroscopic sound 
velocity; Hi) the broadening of the Brillouin peaks, T(Q), 
follows a power law, r(Q) = DQ a , with a « 2 within 
the statistical uncertainties; iv) the value of D does not 
depend significantly on temperature, indicating that the 
broadening (i.e. the sound attenuation) in the high fre- 
quency region does not have a dynamic origin, but rather 
that it is due to the topological disorder Q. 

From a theoretical point of view, two main theories 
have been developed for the study of supercooled and 
glassy systems, one based on first principle computations 
of the equilibrium thermodynamics H and a second one, 
the Mode Coupling Theory (MCT) |j,[7), dynamical in 
nature. Very recently they have been complemented with 
the interpretation of the thermodynamics and of the slow 
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structural dynamics in terms of the topological properties 
of the underlying Potential Energy Surface (PES) [p|-|To[. 
Such interpretation is based on the concept of inherent 
structures, introduced many years ago by Stillinger and 
Weber p]]: upon cooling, the system populates basins 
of the PES whose (local) minima have depth increasing 
on lowering the temperature. Moreover, in this frame- 
work, the dynamics of the systems has been thought to 
be decomposed in a "fast" -high frequency- vibrational 
dynamics describing the exploration of a specific basin 
and a "slow" -diffusive- one, associated to the explo- 
ration of different basins. In these studies a key role has 
been played by computer simulations that, considering 
the nowadays technology, permit a significant sampling 
of the PES also at low temperature and, in principle, the 
study of the dynamics of the model systems at every time 
scale. 

Among other, one not yet fully explained issue in the 
high frequency collective dynamics of supercooled liquids 
and glasses is the recent findings that the temperature 
dependence of the excitations frequency at fixed Q value 
shows a slope discontinuity at a temperature, T Xl that is 
larger than the calorimetric glass transition temperature 
T g [jl2|,[l3). This observation seems to connect to each 
other the shape of the PES basins (which determines the 
vibrational frequency and hence also the excitations fre- 
quency at a given Q's) with the depth of the minima 
explored at that specific temperature. It is our aim to 
further investigate this issue. 

In this paper we show results, obtained by means of 
molecular dynamics simulations, concerning the high fre- 
quency dynamics of a realistic flexible model ]lj-|l6| for 
the fragile glass former orthoterphenyl (OTP). Our aim 
is then twofold. First we want to study the capability of 
the utilized molecular model to reproduce the high fre- 
quency dynamics of OTP by comparing the calculated 
dynamics structure factor with the analogous experimen- 
tal quantity. Once the model has been validate, our aim 
is i) to extend the experimental data to other Q's and 
other thermodynamic points and ii) to provide an inter- 
pretation of the experimental finding of the existence of 
a slope discontinuity in the temperature dependence of 
the sound velocity jl^,|l3|. 

The paper is organized as follows: in Sec. || we briefly 
describe the model and we define the dynamic structure 
factor S(Q,u>) that wi ll b e compared with the experi- 
mental results. In Sec. 
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we recall some of the theoret- 
ical model used in the analysis of the dynamic structure 
factor in liquids; we compare our results with the exper- 
imental sets of data Jl^,[l3| in a large temperature and 
momentum range finding a very good agreement; we dis- 
cuss the connection between our findings and the PES 
approach to the dynamics. Finally in Sec. ^ we discuss 
the results obtained and we draw some conclusions. 



II. COMPUTATIONAL DETAILS 

In our model jl4| the OTP molecule is constituted by 
three rigid hexagons (phenyl rings) of side L a = 0.139 
nm. Two adjacent vertices of the parent (central) ring 
are bonded to one vertex of the two side rings by bonds 
whose length, at equilibrium, is = 0.15 nm. In this 
scheme, each vertex of the hexagons is thought to be 
occupied by a fictious atom of mass Mch = 13 a.m.u. 
representing a carbon- hydrogen pair (C-H). In the iso- 
lated molecule, at equilibrium, the two lateral rings lie in 
planes that form an angle of about 54° with respect to the 
central ring's plane. The three rings of a given molecule 
interact among themselves by an intra-molecular inter- 
action potential; such potential is chosen in such a way 
to give the correct relative equilibrium positions for the 
three phenyl rings, to preserve the molecule from "disso- 
ciation", and to represent at best the isolated molecule 
vibrational spectrum. The interaction among the rings 
pertaining to different molecules is described by a site- 
site pairwise additive Lennard- Jones 6-12 potential, each 
site corresponding to one of the six hexagons vertices. 
The details of the intra-molecular and inter-molecular 
interaction potentials, together with the values of the 
involved constants, are reported in Ref. (m). Previous 
studies of the temperature dependence of the self diffu- 
sion coefficient jl4| and of the structural (a) relaxation 
times Jl4|,[l5| indicate that the introduced model is capa- 
ble to quantitative reproduce the dynamical behavior of 
the real system, but the actual simulated thermodynamic 
temperature has to be shifted by «20 K upward. In the 
following, as is our aim to compare the simulation results 
with the experiments, the reported MD temperatures are 
always shifted by such an amount. 

We have studied a microcanonical (constant energy) 
system composed by 108 molecules (324 rings, 1944 
Lennard Jones interaction sites) enclosed in a constant 
volume cubic box with periodic boundary conditions. To 
integrate the equations of motion we have treated each 
ring as a separate rigid body, identified by the position of 
its center of mass Ri and by its orientation expressed in 
terms of quaternions cji p7[] . The standard Verlet leap- 
frog algorithm |17| has been used to integrate the trans- 
lational motion while, for the orientational part, a refined 
algorithm has been used [Q. The integration time-step 
is 5t=2 fs, which gives rise to an overall energy conserva- 
tion better than 0.01 % of the kinetic energy. 

We studied a wide temperature range spanning the liq- 
uid phase and reaching the glass region. In the different 
temperature runs the size of the cubic box has been cho- 
sen in order to keep the system at the coexistence curve 
experimental density |h]J2(J . At each temperature, after 
an equilibration run (15 ns long) we start the calcula- 
tion of the molecular dynamics trajectory. It is worth 
to note that the system in the low temperature side of 



2 



the supercooled region and in the glassy phase in not 
equilibrated, as the diverging structural relaxation times 
do not allow simulation long enough to reach the equi- 
librium condition. Nevertheless, as we are interested to 
the high frequency properties, i.e. to those properties 
deriving from the intra-basin dynamics of the system, 
we think that it is not really important that the inter- 
basin dynamics was really equilibrated. Actually, also 
in the real experiments, whose outcome we are going to 
compare with the simulation, one investigates the glassy 
phase in similar non-equilibrium situation: the relaxation 
times in the glassy phase are very long, even in compar- 
ison with the experimental measuring times. During the 
MD evolution, lasting for £m=640 ps, the configurations 
of the system are stored for subsequent analysis every 
20 (50) time step for the runs at T > (<) 280 K. All 
the calculations have been performed on a cluster of four 
a-CPU with a frequency of 500 MHz; every nanosecond 
of simulated dynamics needed approximately 24 hours of 
CPU-time. 

From the stored configuration we calculate the dy- 
namic structure factor as measured in an Inelastic X-ray 
Scattering (IXS) experiment @|||. The IXS technique 
measures the electron charge density fluctuation correla- 
tion function. As we consider the phenyl ring as a rigid 
body -i. e. we are not interested in the intra-ring vi- 
brational dynamics that takes place at frequencies much 
higher than those investigated here- we can consider as 
fixed the phenyl ring charge distribution. Within this 
approximation, the scattering center can be considered 
as the ring center of mass, and the effect of the spatial 
distribution of the electron charge is summarized in a 
Q dependent phenyl ring form factor not affecting the 
frequency shape of S{Q 1 oj). Moreover, as we are not 
comparing neither the calculated and measured absolute 
scattering intensities nor their Q-dependence, we neglect 
hereafter the presence of such a form factor. Therefore, 
the appropriate dynamic structure factor plfl to be com- 
pared with the IXS experimental results is given by the 
power spectrum of the phenyl ring center number density 
fluctuation 5pq(t): 
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In this equation tja is the observation time for the vari- 
able Spgit). The dynamic structure factors have been 
evaluated at the Q values allowed by periodic boundary 
condition of the simulation box: Q = 2n/L (n,m,l), 
with n, to, I integers and has been sampled at circular 
frequency between and 22 ps" 1 with a step of «0.075 
ps" 1 . The smallest accessible Q value is ssl.8 nm" 1 . 




FIG. 1. Selected examples of the Q dependence of the 
S(Q,uj) calculated from the MD runs (open circles) at T —50 
K and at the indicated Q values. The data are shown in linear 
and logarithmic scale, in the main figures and in the insets 
respectively. The full lines represent the fit of the data to 
Eq. (Q), while the dashed (dot-dashed) lines is the individual 
elastic (inelastic) contributions to Eq. (Q). 

To increase the statistics, the S(Q,uj) at different Q 
values have been binned in channels 0.2 nm" 1 wide, a 
value comparable to the experimental Q resolution. The 
frequency resolution, dictated by the time extension of 
the runs, is Acj=0.15 (0.75) ps" 1 for T > (<) 280 K. 
The choice of different frequency resolutions below and 
above 280 K is due to the fact that below this temper- 
ature the width of the central line (proportional to the 
inverse of the structural relaxation time r a ) becomes so 
small that is no longer measurable. It is therefore prefer- 
able to relax the energy resolution in order to increase 
the statistics of the data. Above 280 K, the width of the 
central line is measurable, but a "good" energy resolu- 
tion is needed; this procedure produces, obviously, data 
with much poorer statistics. 



III. RESULTS 



Selected S(Q,u) ||] at T=50 K and at Q=1.8, 2.6, 
3.2 and 3.7 nm" 1 are reported in Fig. |l| (open circles). 
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FIG. 2. Selected examples of the T dependence of the 
S(Q,lu) at the fixed Q value of 2.6 nm -1 . Similarly to Fig. |l| 
the data are shown in linear and logarithmic scale in the 
main figures and in the insets respectively. The full lines 
represent the fit of the data to Eq. (0), while the dashed 
(dot-dashed) lines is the individual elastic (inelastic) contri- 
butions to Eq. (Q) . 

In all the figures, to be consistent with the experimental 
data, the S(Q, lu) have been reported as a function of the 
energy (E) measured in meV rather than as a function 
of the circular frequency to. For completeness we recall 
that w[ps _1 ]«1.61 _E[mev]. As the inelastic features ap- 
pear as weak shoulder of an intense central peak, for each 
Q values the S{Q 1 lu) are shown both in linear (main fig- 
ures) and in logarithmic (insets) intensity scale. We are 
interested in extracting the relevant parameters from the 
calculated S(Q,uj). 

A formally exact way to treat these functions goes 
through the description of the density fluctuation cor- 
relation function F(Q,t) = (8pQ(t)5p*Q(Q)), i. e. the fre- 
quency Fourier transform of S(Q,u>) in term of a gener- 
alized Langevin equation [ p3[ : 

F(Q 1 t)+uj 2 F(Q,t)+ f m(Q,t-t')F(Q,t')dt = (3) 



where Lu 2 =kBTQ 2 /M S(Q), m(Q, t) is the "memory func- 
tion" and S(Q) is the static structure factor. This equa- 
tion is "exact", but all the difficulties associated to the 
calculation of the S(Q,lu) have been transferred to the 



determination of m(Q, t). The advantage of this formula- 
tion stays in the fact that -independently from the choice 
of m(Q,t)- the first two sum rules for S(Q,oj) are auto- 
matically satisfied: 

J ' duj S(Q,u)=F(Q,t = Q) = S(Q) (4) 
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duj lu 2 S(Q, lu) = -F{Q, t = 0) 
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By Fourier transform of Eq. (0), it is easy to show that: 



S(Q,cu) = 



S(Q) Lu 2 m'(Q,Lu) 



um"(Q,uj)Y + [um'(Q,u)Y 



(5) 



where m'(Q, lu) and m"(Q, lu) are the real and imaginary 
part of the time Fourier transform of the memory func- 
tion. In the LUT a >>1 limit, a limit that is valid in all the 
investigated temperature range p0| ], the memory func- 
tion can be considered as the sum of two contribution: a 
constant -which reflects the a-process frozen on the time 
scale of the sound waves- plus a very fast decay at short 
time. The latter contribution to the memory function 
-often referred to as "microscopic" or "instantaneous" - 
is usually represented as a delta-function. Introducing 
the two constants representing the area of the "instan- 
taneous" process and the long time limit of the memory 
function, 2T(Q) and A 2 t (Q), the memory function is ap- 
proximated by 

m(Q,t) = 2T(Q)S(t) + A 2 a (Q) (6) 

and therefore, using Eq. (pi), the S(Q,lu) reduces to: 



S(Q,u) = S{Q) x 
f Q S(Lu) + (1 - f Q ) 



n 2 (Q)T(Q) 



^ 2 -ft 2 (Q)) 2 +^ 2 r 2 (Q) 



(7) 



where fi(Q) = v/A 2 (Q)-^ 2 and f Q = 1 - uu 2 jn\Q). 
This expression is the sum of an elastic line (the frozen a 
process) and of an inelastic feature which is formally iden- 
tical to a Damped Harmonic Oscillator (DHO) function; 
the central line accounts for a fraction Jq -the Debye- 
waller or non-ergodicity factor- of the total intensity. 

Recently it has been shown both via MD p3] and via 
IXS [^5| that the "microscopic" contribution to the mem- 
ory function cannot be represented by a delta-function. 
There are, indeed, clear indications that this "micro- 
scopic" part is responsible for a positive dispersion of 
the sound velocity in " harmonic" model glasses [|4j , and 
is responsible for the majority of the positive dispersion 
of the sound velocity in liquid lithium 25 1 . The simplest 
generalization of the delta function appearing in Eq. (|^) 
considers a simple Debye-like " microscopic" contribution 
to the memory function: 



m(Q,t) = A 2 (Q)e 



+ A£(Q). 
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FIG. 3. Q dependence of the excitations frequency Q(Q) 
at the indicated selected temperatures (open symbols for MD 
data and closed symbols for IXS data) as derived from the fit 
of the S(Q,ui) to Eq. (R). The dashed lines are guide to the 
eye. 
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FIG. 4. Q dependence of the excitations broadening T(Q) 
at the indicated selected temperatures (open symbols for MD 
data and closed symbols for IXS data) as derived from the 
fit of the S(Q,u) to Eq. (§). The full line is a quadratic 
(F(Q) = dQ 2 ) fit to the whole set of data. 



The quantity entering in this equation is the char- 
acteristic time of the "microscopic" contribution to the 
memory function. Obviously, whenever J7t m << 1, 
Eq. (H) produces for S(Q,u>) a DHO-like lineshape with 
fi(Q) = y/Al(Q)+u,* and T(Q) = ^{Q)^. 

In analyzing the IXS spectra of OTP, the author of 
Ref. |Q used the simplified DHO expression for the 
S(Q, uj) without noticing any systematic deviation of the 
experimental data from the fitting lineshape. 
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FIG. 5. T dependence of the excitations broadening T(Q) 
at the indicated selected Q values as derived from the fit of 
the S(Q, cj) to Eq. (|). 

Here, in order to better compare the MD results with 
the IXS one, we will follow the same procedure and ana- 
lyze the calculated S(Q,uj) to the lineshape reported in 
Eq. (0). It is worth to note, however, that the MD data, 
less noisy and less affected by the resolution function than 
the experimental data, actually allows to establish that 
the S(Q, oj) obtained by the use of the memory function 
in Eq. (@) produces a better fit of the data. The fits 
to the MD S{Q 1 ui) using Eq. (0) has been made via the 
minimization of the standard \ 2 function, and the results 
are reported in Fig. |l| as full lines. The individual con- 
tributions (elastic and inelastic) are also shown (dashed 
and dot-dashed lines). The agreement between the MD 
data and the fitting function is satisfactory for all the in- 
vestigated T and for all the investigated Q. An example 
of the T dependence of the S(Q,u>) at Q = 2.6 nm _1 
is reported in Fig. ^, again in linear (main figure) and 
logarithmic (inset) scale and together with the best fit to 
a DHO lineshape. 

The Q-dependences of the DHO parameters ft(Q) 

and T(Q) together with the experimental results of 

Refs. [ 12|Jl^1 at selected temperatures are reported in 
Figs. I and || respectively [^6|. These figures confirms 
the experimental findings (see Fig. 1 in |12| and Fig. 2 
in that: i) The excitations frequencies f2(Q) show 
a clear Q dependence, which approach a linear behavior 
at small Q's and bent down at increasing Q values; ii) 
The slope at small Q of il(Q) shows a marked T depen- 
dence together with its overall shape; iii) The parameters 
r(Q) follow a Q 2 behavior; and iv) The temperature de- 
pendence of r(Q) is much less pronounced than that of 
fi(Q). 

The latter T dependence is better emphasized in Fig. ||, 
where the parameter T(Q) is reported -for the indicated 
selected Q values- at all the investigated temperatures. 
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FIG. 6. T dependence of the excitations frequency Q(Q) 
at Q—2.6 nm -1 (scaled at Q=2.5 nm -1 as explained in |E7l) 
and as derived from the fit of the calculated S(Q, u>) to Eq. (pi) 
(open circles). Also reported are the values of Q(Q) at Q—2.5 
nm -1 determined from the IXS experiment (full squares) Jl2| . 
The full lines are the best fit to the experimental points in the 
low- and high-temperature regions: they cross at the tempera- 
ture T x . The dashed lines indicate the limit of ±<r predictions 
band. 

A decrease (line narrowing) of T(Q) of about 20% is ob- 
served between the glassy phase at 50 K and the liquid 
phase (450 K). This decreasing behavior can be explained 
remembering that, in the f2(Q)r M << 1 limit, the param- 
eter T(Q) is given by: 



r(Q) = A*(Q)r M = (nl(Q) - n 2 (Q))r» 



(9) 



and considering that both il(Q) and Qoa(Q) are linear in 
Q and proportional to the sound velocities appropriate 
to frequencies always much larger than 1 Jr a but respec- 
tively smaller and larger than The decrease of T(Q) 
with increasing T reflects therefore the decrease of the 
sound velocities with increasing T. It is also important 
to point out that Eq. (^) gives also an explanation for the 
observed Q 2 behavior of T(Q), when the reasonable hy- 
potheses of a weak Q dependence of r M and of the sound 
velocities are made. Overall the MD results reported in 
Figs. § - ^| are in agreement with the IXS experimental 
findings and give further support to the picture where 
propagating high frequency sound modes exist in glasses 
and liquids. These excitations are characterized (in the 
small Q limit) by a linear and a quadratic Q dependence 
of the excitations frequency and sound absorption coeffi- 
cient respectively. 

To be more quantitative, we report in Fig. ^| the com- 
parison of the T dependence of the £l(Q) parameter at 
fixed Q value of 2.5 nm -1 |^7| as derived from the DM 
and IXS data. The two set of data are in satisfactory 
agreement. 
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FIG. 7. T dependence of the excitations frequencies Q(Q) 
at the indicated Q values derived from the fit of the calcu- 
lated S(Q,lu) to Eq. (Q) (open symbols). The full lines are 
the best fit to the points in the low -and high- temperature 
regions: they cross at the temperature T X (Q), indicated by 
stars. The dashed lines indicate the limit of ±<r predictions 
band: their lowest and highest T crossing are used to define 
the error bar for T X (Q). The dot-dashed line is a guide to the 
eye and connects the different discontinuity temperature. 



The MD data lie slightly above the experimental ones, 
a feature that we ascribe to the non perfect tuning of 
the molecular potential model. It is important to un- 
derline that, similarly to the experimental data, also the 
MD derived fl(Q = 2.5 nm -1 ) show a discontinuity of 
slope at a temperature around T g . A discontinuity that 
also in the MD data appears to be located at a tem- 
perature slightly above the calorimetric glass transition 
temperature. However, the single set of data at Q = 2.5 
nm -1 is not conclusive. A more complete picture can 
be derived from the data reported in Fig. ^, where the 
T dependence of f2(Q) is reported for four different Q 
values. Here one can clearly see that the temperature 
T X {Q) where is observed the slope discontinuity shows 
a well defined Q dependence, and seems to approach T g 
decreasing Q. This observation, obviously, indicates that 
the shape of the dispersion relation fl(Q) is modified by 
changing T. Therefore it is not possible to scale all the 
f2(Q) vs. Q curves on top of each others by a factor, 
and it is important ot study the Q—>0 limit of these 
dispersion relations. This extrapolation is illustrated in 
Fig. H where, at four selected temperatures, the "appar- 
ent" sound velocity c(Q) = Vl(Q)/Q is reported vs. Q. 
A fit of these apparent velocities to a quadratic function, 
c(Q) = c(Q = 0) + aQ 2 is then performed (dashed lines 
in Fig. 8). The Q — sound velocity, obtained by such 
an extrapolation are reported in Fig. ^, together with the 
similar quantity derived from the IXS data p3[. 
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FIG. 8. The apparent sound velocities c(Q) = Q(Q)/Q are 
reported as a function of Q at the four indicated temperatures. 
The dashed lines are the best fits of c(Q) to a quadratic func- 
tion, c(Q) = c(Q = 0) + aQ 2 , and are used to determine the 
Q—0 limit of the sound velocity. 
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FIG. 10. Q dependence of the discontinuity temperature 
T X (Q) derived from Fig. Q (open circles) and from Fig. |^ (open 
square) . The error bars are determined by the crossing of the 
±cr prediction bands. The dashed line is a guide to eye. The 
orizontal arrow indicate the calorimetric glass transition tem- 
perature. 
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FIG. 9. T dependence of the sound velocity at Q=0 de- 
rived from the fits as in Fig. || (open circles). Also reported 
are the corresponding values determined from the IXS ex- 
periment (full squares) [ jUjj . The full lines are the best fit 
to the MD points in the low- and high-temperatue regions: 
they cross at the temperature T X (Q — 0). The dashed lines 
indicate the limit of ±cr predictions band. 

Again we found a good agreement between simulated and 
experimental data, and again the two sets of data show 
a slope discontinuity at a comparable temperature, that 
now coincide within the statistic uncertainties with T g . 

The overall picture is summarized in Fig. [HJ where the 
discontinuity temperature T X (Q) is reported as a function 
of Q. Despite the large statistical error bars, it is evident 
that T x has a marked Q dependence and, in particular, it 
is clear that it extrapolates to a temperature compatible 
with T g as Q goes to zero. 



This result can be interpreted in the framework of the 
current understanding of the dynamics of supercooled 
liquids. Numerous MD simulations study [j^-fiof have 
recently pointed out the validity of the description of 
the supercooled liquid state in terms of the formalism 
initially introduced by Stillinger and Weber jnj: the 
point that represent the system in the 37V configura- 
tional space moves, in this space, performing a "fast" 
vibrational dynamics around a given minima of the po- 
tential energy hyper-surface and occasionally -with a rate 
dictated by the structural relaxation time r Q - perform 
jumps between basins pertaining to different minima 
This "naive" description acquires more physical meaning 
with the demonstration ]28|-|3l|] that the two dynamical 
regimes are almost decoupled, to the extent that the ag- 
ing processes in glasses can be quantitatively described 
supposing that the two subsystem are equilibrated at 
different temperatures. In this framework, it has been 
demonstrated [ |30pl| ] that, approaching the glass transi- 
tion temperature, the representative point visits minima 
of lower and lower potential energy, and exist a strict 
relation between the temperature and the energy of the 
minima visited. As demonstrated in the case of Binary 
Mixture Lennard Jones (BMLJ) (3^,0 systems, different 
minima have also -on average- different metric properties, 
as for example the curvatures at the minimum point and 
therefore the distribution of vibrational eigenfrequencies. 

Our findings on the Q-dependence of the discontinu- 
ity temperature T x can be interpreted within this frame- 
work. Similarly to Lennard- Jones systems, also in the 
case of OTP we found that the modification of the eigen- 
frequencies (local curvatures) does not follow a simple 
scaling, but there is a deformation of the shape of the 
density of vibrational states in changing the energy of 
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the minima. In particular, upon cooling, i. e. on going 
to lower and lower energy minima, the highest curvatures 
(highest frequency, highest Q) reach their limiting values 
(those pertaining to the "ideal" glassy minima) before 
than the lowest curvatures do. 



Finally, our results indicate a possible experimental way 
to validate or disprove the PES-based interpretation of 
the glass transition phenomenology. 
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IV. CONCLUSIONS 

In this paper we have reported results about Molecular 
Dynamics calculation of the high frequency dynamics on 
the flexible molecule model of orthoterphenyl introduced 
in Ref . JTJ] . The collective dynamics structure factor has 
been determined in a wide exchanged momentum and 
temperature range and has been compared to the its ex- 
perimental measurements performed by means of Inelas- 
tic X-ray Scattering |l2|,|l3|,|o| . The relevant parameters 
that describe the overall spectral shape have been deter- 
mined using the same model already utilized in the anal- 
ysis of the experimental data, i. e. the sum of an elastic 
contribution and a DHO lineshape. The resulting param- 
eters fl(Q) and T(Q) -the position and broadening of the 
inelastic peaks respectively- are in good agreement with 
the corresponding experimentally determined quantities. 
Moreover the present simulation confirms the existence 
in liquids and glasses of propagating sound waves at high 
frequency and shows that: i) The excitation frequencies 
Q(Q) show a clear Q dependence, which approach a lin- 
ear behaviour at small Q's and bent down at increas- 
ing Q values; ii) The slope at small Q of Q(Q) shows 
a marked T dependence together with its overall shape; 
hi) The parameters T(Q) follows a Q 2 behaviour; and iv) 
the temperature dependence of the high frequency sound 
attenuation coefficient, related to T(Q), is only weakly 
dependent on T . The latter result further support the 
findings of a sound attenuation mechanism not related 
to anharmonicity 

More important, we hnd that -similarly to the exper- 
imental findings of Ref. the temperature depen- 
dence of the excitation frequency at constant Q value 
shows a slope discontinuity at a temperature T x close to, 
but definitively larger than, the calorimetric glass tran- 
sition temperature. By extending the Q values where 
this analysis is performed, we found clear evidence of the 
Q-dependence of T x . In particular, we find that T X (Q) 
approaches T g in the Q— >0 limit. The existence of a T de- 
pendence of the S(Q,lu) peaks position can be explained 
by the results of recent MD studies of Lennard-Jones 
systems where the T-dependence of the curvatures of the 
PES at the inherent structures were clearly evidentiated. 
The findings of a Q-dependent transition temperature is 
an indication of the fact that the density of vibrational 
states does not change according to a frequency scaling 
law. Rather it is deformed in such a way that the high- 
est frequency reach the limiting value pertaining to the 
"ideal" glassy minimum before than the lowest frequency. 



We thanks R. Di Leonardo and T. Scopigno for intensive 
discussions during the data analysis. 
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